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ABSTRACT 

An observational approach is presented to constrain the global structure and 
evolution of the intracluster medium based on the PCS AT and ASCA distant 
cluster sample. From statistical analysis of the gas density profile and the con¬ 
nection to the Lx — T relation under the /3-model, the scaled gas profile is found 
to be nearly universal for the outer region and the luminosity evaluated outside 
0.2r5oo is tightly related to the temperature through oc rather than T^. On 
the other hand, a large density scatter exists in the core region and there is clearly 
a deviation from the self-similar scaling for clusters with a small core size. A di¬ 
rect link between the core size and the radiative cooling timescale, fcooi and the 
analysis of X-ray fundamental plane suggest that fcooi is a parameter to control 
the gas structure and the appearance of small cores in regular clusters may be 
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much connected with the thermal evolution. We derive the luminosity-‘ambient 
temperature’ relation (Lx — T'), assuming the universal temperature prohle for 
the clusters with short cooling time and adopting a correction T' = 1.3T, and hnd 
the dispersion around the relation signihcantly decreases in comparison to the 
case of the Lx — T and the slope becomes less steep from 3.0llo;44 to 2.80l^;^t 
-hikev, which is defined as a normalization factor for each cluster, can be re¬ 
garded as constant for a wide range of the cooling time. We further examined 
the Lx — TjS and Lx — T'(3 relations and showed a trend that merging clusters 
segregate from the regular clusters on the planes. A good correlation between the 
cooling time and the X-ray morphology on the Likev — (Looi/fage) plane leads us 
to dehne three phases according to the different level of cooling, and draw a phe¬ 
nomenological picture: after a cluster collapses and fcooi falls below the age of the 
universe, the core cools radiatively with quasi-hydrostatic balancing in the grav¬ 
itational potential, and the central density gradually becomes higher to evolve 
from an outer-core-dominant cluster, which marginally follows the self-similarity, 
to inner-core-dominant cluster. 

Subject headings: galaxies: clusters: general - X-Rays: galaxies: clusters 


1. INTRODUCTION 

The X-ray luminosity-temperature (Lx — T) relation of galaxy clusters is one of the 
most fundamental parameter correlations, established from previous X-ray observations (e.g. 
Edge et ah 1990). Since the X-ray luminosity reflects temperature and density prohles of 
hot intracluster medium (ICM), the Lx — T relation should contain information on physical 
status and evolution of the ICM. 

Observationally, the correlation is well approximated with a power-law function: Lx oc 
T", with a ~ 3 (e.g. Edge et al. 1990; David et ah 1993; Markevitch 1998; Arnaud & 
Evrard 1999). In addition, it shows a signihcant scatter around the mean power-law relation 
and little redshift evolntion. On the other hand, the self-similar model (e.g. Kaiser 1986) 
predicts o; = 2. Thus the inconsistency between the observations and the simple theoretical 
model has been debated for many years and various possibilities including non-gravitational 
heating (e.g. Evrard & Henry 1991; Cavaliere, Menci, & Tozzi 1997) and dependence of gas 
mass or gas-mass fraction on the temperatnre have been proposed (e.g. David et al. 1993; 
Nenmann & Arnand 2001). Recently, hydrodynamical simulations have indicated that the 
effect of radiative cooling plays an important role to reprodnce the observed Lx — T relation 
(e.g. Mnanwong et al. 2002). They also snggested a fnrther reqnirement of signihcant non- 
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gravitational heating mechanism so as to account for the observed gas-mass fraction, however, 
the physical origin of the additional heating is yet to be understood. 

Under the virial theorem and the isothermal /5-model (Cavaliere & Fusco-Femiano 1976), 
an ICM temperature Tgas is proportional to a virial temperature, and differs by a 
factor of /5, Tvir ~ T^asP- However, since Tgas is usually an emission-weighted temperature 
measured from X-ray spectroscopy and strongly reflects the temperature of the cluster core 
region, a temperature decrease due to radiative cooling may have much influence on the 
Tgas measurement. The ASCA spectroscopy of central region of clusters particularly with 
cD galaxy provided important results that cool emission from the core is systematically 
less luminous than that estimated in the past mainly on the basis of imaging observations 
(e.g. Makishima et al. 2001). After the advent of the XMM-Newton and Chandra X-ray 
satellites, more detailed temperature profiles have been measured for nearby ‘cooling-flow’ 
clusters, which typically show central temperature decrement by a factor of ~ 2 (e.g. Kaastra 
et al. 2004). These facts gave a piece of evidence against the standard cooling-flow model 
(Fabian 1994), which triggered explorations of a variety of new scenarios for gas heating: 
heat conduction from outer hot layers (e.g. Narayan & Medvedev 2001), AGN heating (e.g. 
Bohringer et al. 2002), magnetic reconnections (Makishima 1997), Tsunami (Fujita et al. 
2004) etc. Recently Masai & Kitayama (2004) proposed a quasi-hydrostatic model for cluster 
gas under radiative cooling, which predicts a characteristic temperature prohle with an 
asymptotic temperature for the central region being ~ 1/3 of the non-cooling outer region. 
Their calculation agrees well with that observed in the cooling-flow clusters. If the ‘universal’ 
temperature prohle emerges in all or a fraction of clusters with short cooling timescales, the 
effect of temperature drop should be properly taken into account in the discussion of the 
Lx — T relation. 

Another point is due to a connection of the Lx — T relation to the cluster core sizes. 
Ota & Mitsuda (2002) found based on the X-ray data analysis of 79 distant clusters with 
ROSAT and ASCA that the distribution of the cluster core radius shows two distinct peaks 
at ~ 50 kpc and ~ 200 kpc for Hq = 70. When dividing the sample into two subgroups 
corresponding to the two peaks, they show a signihcant difference in the normalization factor 
of the Lx — T relation. Some possibilities have been discussed to understand the origin of the 
double-peaked distribution (Ota & Mitsuda 2002, 2004), however, it is yet to be clarihed. 
As long as one rely on the hydrostatic equilibrium and the isothermal /5-model, their result 
may indicate that underlying dark matter distribution is likely to have two preferable scales 
of 50 kpc and 200 kpc. On the other hand, it is also found that the core radius is tightly 
related to the radiative cooling timescale (e.g. Ota & Mitsuda 2004; Akahori & Masai 2005). 
Thus the core structure might largely be a result of the thermal evolution of the ICM. 
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In this paper, using a large number of ROSAT and ASCA clusters compiled by Ota 
(2001); Ota & Mitsuda (2004), we investigate correlations between the fundamental cluster 
parameters. Since ASCA has a high sensitivity to measure the X-ray spectrum in the wide 
energy band while ROSAT is good at imaging in the soft X-ray band, the two observatories 
are an excellent combination to study properties of the intracluster medium. At present the 
XMM-Newton and the Chandra satellites are in orbit and generate much cluster data with 
higher sensitivities. However, the data set used in the present paper will be one of the best 
existing to construct the largest sample of distant clusters and study global X-ray structures. 
While Ota & Mitsuda (2002) mainly focused on the discovery of the typical two core sizes 
based on the analysis of 79 distant clusters with ROSAT and ASCA^ Ota & Mitsuda (2004) 
described the uniform analysis, and showed the X-ray parameters for the individual clusters 
and the results on scaling relations. The present paper addresses the Lx — T relation in 
more detail, which is organized as below. In § 2, the sample and parameter estimation are 
described. In § 3, we show gas density prohles of the sample and how they are related to a 
scatter about the Lx — T relation. Since redshift evolution is not clearly seen in the present 
distant sample and our analysis is based on a conventional /5-model, a study of redshift 
evolution in the Lx — T relation is beyond the scope of the present paper. In § 4 and § 5, 
we mainly focus on the Lx — T relation and show its connection with the cluster core radius 
and the cooling time. In § 6 we also derive a Lx — Tj3 relation and show a trend that 
possible merging clusters segregate from regular clusters on the plane. Finally in § 7, major 
results are summarized and the evolution of the ICM are discussed particularly in the light 
of radiative cooling. 

We use Hm = 0.3, Ha = 0.7 and h-jo = i7o/(70 kms“^Mpc“^) = 1. The quoted 
parameter errors are the 90% confidence range throughout the paper unless otherwise noted. 
The la error bars are plotted in all figures. 


2. SAMPLE AND PARAMETER ESTIMATIONS 
2.1. The Sample 

The sample is selected from an X-ray catalogue of the ROSAT and ASCA distant clusters 
presented in Ota & Mitsuda (2004). Though the original catalogue includes 79 clusters in 
the redshift range of 0.1 < 2 : < 0.82, the hnal sample used in this work comprises 69 clusters 
in 0.10 < 2 < 0.56 in order to reduce possible selection effect and measurement uncertainties 
which are more serious at higher redshifts, as detailed below. 

Because the sample was collected from the archival data of the pointing observations, in 
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order to reduce the selection effect that is not negligible at higher redshift, we exclude 
5 high-redshift clusters with large measurement uncertainties from the Ota & Mitsuda 
(2004) sample: they are #75 MS2053.7-0049, #76 3C220.1, #77 MS1137.5+6625, #78 
RXJ1716.6+6708, #79 MS1054.5-0321 (see Ota & Mitsuda (2004) for details). We ex¬ 
clude #71 3C295 and #76 3C220.1 due to possible contaminations of X-ray emission from 
the central radio sources as reported by Chandra observations (see Allen et ah (2001a) for 
3C295 and Worrall et ah (2001) for 3C220.1), and #56 CL0500-24, #61 CL2244-02, and 
#66 CL0024+17 due to uncertainties in the spectroscopy caused by contaminations from 
point sources. We did not use #43 A1758 and #79 MS1054.5-0321 because their /5-model 
parameters were not well constrained. 


2.2. Brief Summary of X-ray Data Analyses 

The data analyses consisted of two major steps: 1) the spatial analysis with the ROSAT 
HRI and 2) the spectral analysis with the ASCA GIS and SIS. 

For 1), the radially-averaged X-ray surface brightness distribution was htted with the 
isothermal /5-model to determine the slope parameter (3, the core radius Vc and the central 
surface brightness. The model gave an reasonably good £t to the observed surface brightness 
of many clusters: it resulted in a statistically acceptable fit for about 2/3 of the sample. For 
some clusters, however, it has left residual emission in the central region and a significant 
improvement was found when including an additional /5-model component particularly for 
9 clusters (termed as “double-/5” clusters). We have further found that adoption of the 
single /5-model for such double-/5 clusters will return parameters consistent with one of the 
two /5-model components, as shown in Fig. 9 of Ota & Mitsuda (2004), which means the 
single /5-model tends to pick up a component that dominates the emission. Thus the present 
work basically stands on the single /5-model for the reason that we are mainly interested 
in a systematic study on the global structure in a large number of distant clusters. Note 
that possible systematic errors of physical parameters estimated under the /5-model will be 
discussed in the next subsection. Furthermore, the X-ray morphology was measured with 
using ‘centroid variations’ of the cluster image and the sample was classified into “regular” 
and “irregular” clusters (see §3.2 of Ota & Mitsuda (2004) for details). 

For 2), the emission-weighted temperature and bolometric luminosity were measured 
from the ASCA spectroscopic data under the Raymond-Smith model (Raymond & Smith 
1977), utilizing the XSPEC software version 9.0. kT is the emission-weighted temperature 
determined from the simultaneous fitting to the ASCA GIS and SIS spectra extracted from 
r < 6' and r < 3' circular regions, respectively, and the bolometric luminosity. Lx is calcu- 
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lated within a virial radius rsoo- Here rsoo is defined as the radius within which the average 
matter density is equal to = 500 times the critical density of the Universe at the cluster 
redshift. 

In what follows, we use the X-ray parameters derived under the single /5-model; the 
slope parameter (3, the core radius rc, the central electron density Ueo, and the radiative 
cooling time fcooi as well as spectroscopically measured T and Lx- All quantities are listed 
in Tables 2, 4, and 5 of Ota & Mitsuda (2004). 


2.3. Parameter Estimation under the /5-model and systematic errors 

To perform a reliable estimation of the /5-model parameters, we have carefully checked 
the effect of the ROSAT/WRl PSF on the spatial analysis. As detailed in § 3.4 of Ota 
& Mitsuda (2004), we confirmed that the effect is negligibly small and the parameters are 
well reproduced as long as is larger than 5". We have also checked the consistency with 
the previously published values from the Chandra observations in regard to the /5-model 
parameters (/5 and Tc) and the temperature, and found a good agreement (see § 3.4 and 
§ 4.2 of Ota & Mitsuda (2004)). Due to the large ASCA PSF, it is not easy to measure the 
temperature profile accurately for the distant clusters. We then use the spectroscopically 
determined mean temperature in the parameter estimation. We will further evaluate possible 
systematic errors in estimating physical parameters under the /5-model, Lx and tcooi, which 
are important quantities in the present study. 

• Lx 

To estimate a systematic error of Lx, we first calculate the bolometric luminosity 
within r < 6' in two different ways: 1) convert the 2-10 keV luminosity measured 
with the GIS spectrum into the bolometric luminosity using the emissivity of the 
Raymond-Smith plasma model and the ASCA temperature, and 2) integrate the /5- 
model surface brightness distribution determined with the ROSAT HRI within the 
same integration area and multiply it by the band correction factor. We found that 
two estimations are consistent within about 10-20% in all three different temperature 
ranges of 2-5, 5-7, > 7 keV or four different redshift ranges, 0.1-0.2, 0.2-0.3, 0.3-0.4, 
> 0.4. This also assures that the current assumption of the isothermal /5-model gives 
reasonably a good approximation to the real clusters in estimating Lx- Then we finally 
determine the bolometric luminosity within rsoo. Lx from 1) the bolometric luminosity 
that was spectroscopically measured within 6' under the Raymond-Smith model and 
corrected for the emission from outside 6' by assuming the /5-model surface brightness 
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distribution. In Fig. la we show a comparison of the bolometric luminosities measured 
with GIS and HRI. 


In the temperature range concerned, the cooling rate of the gas is dominated by free-free 
emission and then approximately proportional to y/T. Therefore, with the cooling rate 
in the form q^y/TnlQ, we dehne the cooling time for a cluster as tcooi = ^ky/T/q^Ueo, 
where Ueo is the gas density at the cluster center. On the other hand, recent higher- 
resolution observations of clusters with a cooling core pointed that the conventional 
/3-model does not sufficiently fit overall prohles (e.g. Lewis et ah 2003; Hicks et al. 
2002). In addition, if the central temperature drop by a factor of ~ 2 exists like what 
observed in many cooling-flow clusters, it might cause an overestimation of fcooi- We 
then make a comparison with recent results from Bauer et al. (2005) as they analyzed 
Chandra data of X-ray luminous clusters at 0.15 < 2 : < 0.4 and derived the cooling 
time prohles using spatially resolved spectral information. From Fig. lb, we hnd that 
there is an agreement between the two analyses: for 11 out of 21 distant clusters tcooi 
is consistent with the Chandra results within the measurement errors. For the rest of 
them, it is systematically larger (smaller) for fcooi ^ 2 Gyr (fcooi ^ 10 Gyr). As the 
clusters with tcooi ^ 10 Gyr are mostly irregular ones, the underestimation might be 
due to their complex internal structure and different dehnitions of the centroid position. 
From temperature maps derived by Bauer et al. (2005), we found that clusters with 
systematically smaller tcooi tend to have hot spots at the centroids or gas hotter than 
the surrounding regions. This would then bias the temperature toward lower value 
(as we use the emission-weighted temperature instead of the spatially-resolved central 
temperature) and then higher central density (to account for the observed emission 
measure), resulting in shorter tcooi- For clusters with tcooi ^ 2 Gyr, the overestimation 
may be attributed to possible temperature decrease at the center: assuming that the 
central temperature isl/2-l/3of ambient temperature, T', and the emission-weighted 
temperature is related through T ~ T'/1.3 (Eq. 21 in § 5.2), tcooi is possibly overesti¬ 
mated by a factor of ~ 1.3 — 1.8. Thus this may reconcile the factor of ~ 2 difference 
seen in Fig. lb. However, because the above estimation depends on the assumed tem¬ 
perature prohle, we use tcooi obtained from the /3-model analysis of the ROSAT and 
ASCA data. Note that the two analyses compared in Fig. lb are consistent within at 
most a factor of 2 and our discussion is not so much affected by the absolute value 
of tcooi as this paper aims at a study of a global trend of the IGM evolution, which is 
suggested to proceed along the “principal axis” (§ 4.4) over two orders of magnitudes 
of cooling time. 



3. Lx - T RELATION AND GAS DENSITY PROFILES 


3.1. Core Radius Distribution and rsoo — Tc Relation 

Since the Lx — T relation is sensitive to the structure of the cluster core region, we 
hrst examine the core radius. As noted in Ota & Mitsuda (2002), together with 45 nearby 
clusters analyzed in Mohr, Mathiesen & Evrard (1999), the Vc distribution of 121 clusters 
exhibits high concentrations around 50 h^Q^kpc and 200 /ifg^kpc . We show the core radius 
distribution of the present sample in Fig. 2b. Though the has a large cluster-to-cluster 
dispersion, we hud that about 60% of clusters are included in narrow ranges corresponding 
to the two peaks. Thus the large dispersion is partly due to the double peak structure of 
the Tc distribution. 

Under the self-similar model, the internal structure of the gas should be scaled by 
the virial radius, and then t^qq/tc should be constant for all clusters. However, as is seen 
from Fig. 2a, Vc does not simply scaled by rsoo. As mentioned in Ota & Mitsuda (2004), 
we should be careful about the results of the rsoo — Tc relation because it depends on the 
choice of the range. When we £t all the sample at the same time, we obtain a weak 
correlation, rsoo oc (y^/d.o.f. = 953.2/67). Since t^qq/tc 3> 1 in general, rsoo is 

mostly determined by the temperature. Actually Tc and T do not show a strong correlation. 

0 . 37 +^-^^ 

If we divide the sample into two Vc groups relative to Vc = 100 kpc, rsoo oc rc “ for 34 
clusters with Vc > 100 kpc (y^/d.o.f. = 107.7/32) while no meaningful correlation was found 
for 35 clusters with Tc < 100 kpc (the correlation coefficient is 0.21). We further restrict 

the large core group to a narrow range, for example 100 < Tc < 200 kpc, rsoo oc Vc 
(x^/d.o.f. = 21.0/16). Thus we consider that the large core clusters may marginally satisfy 
the self-similar condition (see also Akahori & Masai (2005)). On the other hand, since 
the small core group clearly shows a deviation from the rsoo ~ U relation, the formation 
of the small cores is suggested to be determined from some physical process other than 
the self-similar collapse. Furthermore, for 6 clusters with a very large core radius of Tc > 
400 kpc, which forms the third peak around 500 kpc in the r^ histogram (Fig. 2b), there 
is no meaningful rsoo ~ correlation (the correlation coefficient is 0.21). Thus they do not 
satisfy the self-similarity. This deviation may be attributed to the past or on-going merging 
events, which we will discuss in more detail in § 6.1. 
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3.2. Scaled Gas Density Profiles 

To investigate the global gas density structure, we plot the best-£t /d-model profiles for 
69 clusters in Fig. 3, where the radius was normalized by rsoo- A significant scatter of the 
density exists within O.lrsoo region (typically O.lrsoo 100 kpc) and the smaller core 
clusters show a systematically higher density in comparison to the large core clusters. On 
the other hand, the profiles are fairly universal outside that radius. We consider that this 
is consistent with the similarity found in outer regions of scaled prohles for nearby clusters 
(Neumann & Arnaud 1999). 

In the next two subsections we will treat the core region and the envelope separately 
in order to directly show how the density profile is related to the scatter around the Lx — 
T relation. We found that the rsoo — Tc relation and the density profile do not show a 
clear evolution against the observed redshift. However it is important to test whether the 
Lx — T relation would evolve with redshift or it is more closely related with other physical 
parameters. Then in § 3.5 we will investigate the redshift dependence of the Lx — T relation 
since the present sample covers the redshift range of 0.1 < < 0.56. 


3.3. Lx — T outside 0.2r5oo 

We calculate the luminosities outside a radius of 0.2r5oo and plot them against the 
temperature in Fig. 4b. The radius was chosen so as to be a typical Tc of the large core 
clusters, ~ 200 kpc and thus the emission from the core region can be removed from the 
overall emission. As a result, we obtained 

Lx(> 0.2r5oo) = 4.95^1^3 x 10^^(A;T)2-®®-o:“. (1) 

We found that there is no significant difference in the relation between the two Tc groups (see 
also Table 1) and the data scatter around the mean relation is surprisingly small. Despite 
the similarity of the outer profiles shown in §3.2, however, the measured power-law index is 
yet close to 3 rather than the self-similar expectation of 2. This is also true for 20 clusters 
with 100 < Tc < 200 kpc, which marginally satisfy the self-similar relation as noted in § 3.1; 
the power-law index is 3.13 ^q;5®. 

Note that the Lx — T relation is highly sensitive to the way the central emission is 
treated. Allen & Fabian (1998) pointed out that the slope of the Lx — T relation would be 
flattened from to once corrected the central emission based on the standard cooling 
flow model. We find that the slope is still steeper than the prediction of the self-similar 
model by simply excluding the central component in a model-independent manner. 
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3.4. Lx — T inside 0.2r5oo 

In Fig. 4a, we show the Lx — T relation calculated for r < 0.2r5oo. The power-law htting 
resulted in 

Lx{< 0.2r5oo) = 3.02l^:'8f x . (2) 

There is a signihcant data scatter in the plot and the luminosity is systematically higher 
for the small core group for a fixed temperature. If assuming the best-fit slope of 4.37, the 
relations for the small and large core groups are obtained as: 

Lx(<0.2r5oo) = 5.82lJ;?^ X for re < 0.1 Mpc, (3) 

Lx(< 0.2r5oo) = 8.32+Jj® x for > 0.1 Mpc (4) 

(see also Table 1). This signihcant difference in the normalization factor can be attributed 
to a negative correlation between the central density and the core radius as mentioned in 
§3.2. 

We then see from Figs. 3 and 4 that the density scatter inside 0.2r5oo is indeed a source 
of scatter seen in the Lx — T relation. In comparison to the result for r > 0.2r5oo, the large 
variety of the luminosity for r < 0 . 2 r 5 oo and the signihcant deviation from the self-similarity 
for the small core clusters may imply that the dispersion around the Lx — T is closely related 
to the evolution of the ICM in the central region after the collapse. 


3.5. Redshift Dependence of Lx — T and Comparison to Previous Results 

We investigate the redshift dependence of the Lx — T relation by evaluating the relations 
for the following three subsets of data: 18 clusters in 0.1 < 2 : < 0.2, 27 clusters in 0.2 < 2 : < 
0.3, and 24 clusters in 0.3 < < 0.56 (only 7 out of 24 have > 0.4) and compare with the 

previous results on high-redshift {z > 0.4) and low-redshift {z < 0.1) clusters. To constrain 
the redshift dependence in the core and outer region separately, we also derive the Lx — T 
relation for r < 0.2r5oo and r > 0.2r5oo. We did not apply the evolution correction to Lx (i.e. 
multiplying Lx by E{z)~^ since it does not signihcantly change the relation for Ac = 500 as 
shown in Table 9 of Ota & Mitsuda (2004).) 

• Lx — T relation 

The power-law fitting of the Lx(< rsoo) — T relations for three redshift subgroups 
resulted in 


Lx — 


1.23lg|^ X for 0.1 < z < 0.2, 


( 5 ) 
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Lx = 5.13;1®23 X for 0.2 < z < 0.3, (6) 

Lx = 2.95tl%^ X 10^^(A;T)^-®^-i' 86 for 0.3 < < 0.56 (7) 

(Fig. 5a). There is a marginal trend to obtain a steeper slope and a smaller normaliza¬ 
tion factor for the higher redshht clusters. However, because of the large uncertainties 
seen in the slope parameter, if hxed at the best-£t value of the total sample, 3.01, the 
normalization factors overlap with each other within the errors as shown in Eqs. 8-10 
and Fig. 5b, and thus the redshift dependence is not clear in the observed redshift 
range. 


Lx = r.OAtlil X for 0.1 < < 0.2, (8) 

Lx = 7.08l?;^J X 10^\kTf-^^ for 0.2 < < 0.3, (9) 

Lx = 0.02t\fo X l0*‘^{kTf-^^ for 0.3 < < 0.56. (10) 

Ettori et ah (2004) reported a steep slope of 3.72 ± 0.47 (la error) for high redshift 
clusters (0.4 < 2 : < 1.3) based on the Chandra data analysis and suggested a negative 
evolution. Eqs. 6 and 7 agree with their relation within the errors. Thus there might be 
a weak redshift evolution in the Lx—T relation. However, due to the large uncertainties 
associated with the relations for the higher redshifts, it is not statistically signihcant. 

• Lx(< 0.2r5oo) — T and Lx(> 0.2r5oo) — T relations 

We obtain the Lx(< 0.2r5oo) — T relation in the three redshift ranges as follows (see 
also Fig. 5c). 

Lx(<0.2r5oo) = 3.89l?;| X 10^^(A;T)^-3^ for 0.1 < < 0.2 (11) 

Lx(< 0.2r5oo) = 2.39to}l x l0^\kT)^-^'^ for 0.2 < z < 0.3 (12) 

Lx(< 0.2r5oo) = x lO^^kTY'^'^ for 0.3 < < 0.56. (13) 

There is only a weak signature of negative evolution regarding the Lx — T within 
central 0 . 2 r 5 oo region: this change in the normalization factor against redshift is not 
statistically signihcant and smaller than that of two subgroups with different core sizes 
in the above analysis (Fig. 4a). Furthermore, as shown in Fig. 5d, if the central 0.2r5oo 
region is excluded, the normalization factor of the Lx(> 0.2r5oo) —T relation has little 
redshift dependence. 

In order to further compare our results to the previous results, we followed the Marke- 
vitch (1998)’s analysis as much as possible: the X-ray luminosity was computed within 
a fixed radius of 2 Mpc and the possible cool emission inside 100 kpc was corrected by 
multiplying the luminosity for 0.1 < r < 2 Mpc by 1.06 for clusters with short cooling 
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time, logtcooi/^age < —0.5, while their ambient temperature, T' was estimated by using 
Eq. 21. As a result, there is no signihcant difference in the Lx — T among the three 
subsets of data and they are consistent with the local relation (Fig. 5e). Therefore the 
redshift dependence is not clearly seen in the present sample. 

There are recent Chandra and XMM-Newton studies of high-redshift {z > 0.4) clusters 
which support a positive evolution (Vikhlinin et ah 2002; Kotov & Vikhlinin 2005; 
Maughan et ah 2005; Lumb et al. 2004). We should note that the current result does 
not rule out the scaling of Lx(l + z) found by Vikhlinin et al. (2002): we found 
that the normalization factors of the Tx(l + z)~^'^ — T relations are again consistent 
with that of local relation within their errors (Fig. 5f). In his recent review, Voit 
(2005) pointed that there are actually some systematic uncertainties in the evolution 
measurement and it is not precise enough to distinguish the possibilities of positive or 
negative evolution. 

Based on the above analyses, we may consider that the redshift dependence is likely to 
be a minor effect in discussing the source of scatter around the Lx — T relation. Although the 
marginal, negative evolution of the normalization factor seen in the Lx{< 0 . 2 r 5 oo)—T relation 
might imply a progress of gas concentration as redshift, the study of redshift evolution 
requires higher-accuracy determination of parameters and thus is beyond the scope of the 
paper. The scatter seems to link more strongly to the other parameters describing the core 
structure, which we will try to show more directly in § 4 and § 5. 

Besides, Ota & Mitsuda (2004) showed that there is not strong redshift evolution in the 
observed X-ray properties particularly at z < 0.5, including the central gas density, the core 
radius, and temperature. As we will show in § 4.4, the current sample is also found to have 
an X-ray fundamental plane similar to that found in nearby clusters (Fujita & Takahara 
1999). We will then carry out statistical studies regardless of their redshifts in the paper. 


4. Lx-T RELATION AND 
4.1. Observed Lx — T Relation 

In Fig. 6a, we show the Lx — T relation of 69 clusters where we again divide the sample 
into two subgroups according to two different Tc ranges. As pointed out in Ota & Mitsuda 
(2002), the distributions of the two subgroups are systematically shifted in the hgure; the 
luminosity of a small core cluster is on average higher than that of a large core cluster for 
a hxed temperature. If we £t all the 69 clusters at the same time, we obtain a power-law 
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index of 3.01]'^o;44 (90% error). Since we found that there is not signihcant difference in the 
power-law indices between the two subgroups (Table 1), we obtain the best-£t normalization 
factors of the relation to be l.OStoii x 10^^ ergs“^ and 4.521q;34 x 10^^ ergs“^ by fitting them 
separately but fixing the indices at the mean value of 3.01. Thus the difference is significant: 
the cluster with a smaller core radius is in general brighter for a given temperature. The 
fitting results are summarized in Table 1. 


4.2. Lx — T — Tc Relation 

The above observational results show that the dispersion around the mean Lx — T 
relation is connected with the cluster core size. Assuming the best-fit power-law index of 
a = 3.01, we scale the luminosity of each cluster to the luminosity for temperature of 1 keV 
as Likev = Tx(/cT/lkeV)“". In Fig. 6b, we plot Likev as a function of Vc- We found that 
there is a negative correlation approximated as Tikev oc Vc from the fitting. This 

suggests that Lx depends not only on T but on Vc as 

Lx oc (14) 

The result is consistent with that derived by Fabian et al. (1994), although they interpreted 
the dependence on as the cooling-flow rate. Moreover, Eq. 14 is also consistent with the 
X-ray fundamental plane of nearby clusters noted by Fujita & Takahara (1999). In the next 
subsection, we will derive a simple model for the Lx — T relation under the assumption of 
/5-model, incorporating the two characteristic core radii, and compare with the observation. 


4.3. Comparison to a Lx — T — Tc Model 

We obtained the following observational facts so far: 1) Tc distribution shows high 
concentration at 50 and 200 kpc, which result in a large dispersion of Tc, 2) Vc does not show 
strong correlation with the virial radius, and 3) Lx depends on both T and Vc- These results 
suggest that Vc is determined by some physical process in the core region of the cluster, 
independent of the cluster virial mass. We thus modify the assumption of self-similarity and 
consider its consequence on the Lx — T relation. 

In order to evaluate the Lx — T relation quantitatively, we assume the isothermal (3- 
model for the gas density and a constant gas-mass fraction. The details of the calculation of 
this simple ‘^Lx — T — T'c model” will be mentioned in Appendix A. For a typical value of 
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/9 = 0.7, we obtain the model in the form of: 

Lx oc (15) 

This is close to the observed relation (and the equation for the X-ray fundamental plane, 
Eq. 14), however, the index for Tgas is smaller than the observation. In Fig. 7, we show the 
model curves for three different values of Tc- Thus the Lx — T — model based on the 
isothermal (3 model may roughly fit the observed trends of the Lx — T relation, though not 
sufficient to reproduce the data quantitatively. 


4.4. X-ray Fundamental Plane Analysis 


To look further for an essential parameter to understand the dispersion around the 
Lx — T relation, we investigate the X-ray fundamental plane in the (logrieo, logT, logrc) 
space for the distant sample. We utilize a method similar to Fujita & Takahara (1999) but 
carry out a £t to take into account the statistical errors. As a result, we obtained three 
orthogonal parameters as follows: 


X 

Y 

Z 


oc 

oc 

oc 


^0.44 0.65^-0.62 
"'eO ' c ) 

0.45^0.44^0.78 
"'eO ' c ) 

„0.78„-0.62Tn-0.10 
'^eO ' c ^ 


(16) 

(17) 

(18) 


In Fig. 8a, the log X —log Z plane for the distant clusters is shown. We found that the above 
equations agree with those obtained for the nearby sample by Fujita & Takahara (1999) and 
then the > 0.1 and < 0.1 clusters are distributed approximately on the same fundamental 
plane. In Fig. 8b we also show the distribution of nearby clusters projected onto the same 
log X — log Z plane, whose log Z range coincides well with that of distant clusters. 

As y^/d.o.f. = 492.1/66 for the first fit (Eq. 16), however, the constancy of log X i.e. the 
planarity of the cluster distribution in the 3D space is rejected in the strict sense. Besides 
the data scatter, this is likely due to the presence of a “break” of the fundamental plane that 
is seen in Fig. 10b. We will go back to this problem in § 5.2. 

By setting X ~ const., Eq. 18 yields 

Z oc oc (19) 

Thus the principal axis of the fundamental plane, Z is suggested to be closely related to r^- 
Since fcooi oc VT/ueo, Eq. 19 can be rewritten as 


( 20 ) 



15 


Since the Z-axis represents a direction along which the dispersion of the data points becomes 
the largest in the parameter correlations, tcooi is considered to be a key parameter to control 
the ICM evolution. Hence we will focus on the effect of gas cooling in the next section. 


5. Lx - T RELATION AND tcooi 
5.1. Radiative Cooling and Temperature Decrease 

As already mentioned in the previous section, there is a large dispersion of the gas 
density inside O.lrsoo region which is strongly related with Tc (Fig. 3). Ota & Mitsuda 
(2004) derived the rieo — Tc relation to be Ueo oc r~^'^ and noted that the correlation tends to 
become even steeper for small core clusters: Ueo oc for Tc < 100 kpc. For clusters with 
such high density cores, radiative cooling may work considerably. It is also suggestive that 
O.lrsoo is roughly equal to the typical cooling radius of rcooi ~ 100 kpc. 

In Fig. 9a we show the radiative cooling time against the core radius. We conhrmed 
that for all small core clusters, it is shorter than the Hubble time (tn = 13.4 Gyr; Spergel 
et ah (2003)), tcooi < ^h- Even if we use the age of the Universe at the cluster redshift tage 
instead of tn, which ranges 8.0 < fage < 12.2 Gyr for the present sample, that hardly alter 
the result. Furthermore considering the strong dependence on the core radius, tcooi oc for 
Tc < 100 kpc (Ota & Mitsuda 2004), these results seem to indicate a central concentration 
of the gas according to the progress of radiative cooling. This is supportive of an idea of the 
standard cooling-flow model in this regard. However, as seen from Fig. 9b, there is no clear 
difference in the temperature range between the two subgroups of different core sizes. The 
average temperature (and the standard deviation) is 5.4 keV (1.8 keV) for the clusters with 
Tc < 100 kpc and 7.0 keV (1.8 keV) for > 100 kpc, and the difference is only 30%. Thus 
there is not a strong temperature drop as predicted by the standard cooling-flow model. 
Our result conhrms the lack of strongly cooled gas based on the large sample of ROSAT and 
ASCA clusters. 


5.2. Lx — T — tcooi and Lx — T' — tcoo\ Relations 

To investigate the nature of the dispersion around the Lx — T relation in more detail, 
we show Likev as a function of the cooling time normalized by the age of the universe at 
the cluster redshift, (Gooi/Gge) in Fig. 10b. At log(tcooi/tage) > —0.5 the distribution is 
nearly constant within the data scatter. On the other hand, at the shorter cooling time 
of log(fcool/Gge) ^ —0.5, a signihcant deviation from the mean Lx — T relation (which is 
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denoted with the horizontal solid line in the hgure) is clearly visible. We found that the result 
depends little on the redshift, as shown in Fig. 11. Although this was interpreted as central 
excess emission accompanied by the radiative cooling within the framework of the standard 
cooling-flow model, we want to suggest another possibility that we might underestimate the 
temperature for the clusters with short cooling time because of mild temperature decrease 
at the center. 

We then evaluate the emission weighted temperature of the clusters with short cooling 
time utilizing the universal temperature prohle found in many nearby cooling-flow clusters. 
Assuming a projected temperature prohle, T(r) oc 0.40 -|- 0.61[(a;/a;c)^'®/(l -|- {x/xcY'^] [x = 
r/r 25 oo) for relaxed clusters (Allen et al. 2001c) and the /? prohle, S'(r) oc (1-|- 
T = T{r)S{r)2nrdr/ S{r)27irdr . Since the temperature prohle is nearly constant 
outside the typical cooling radius of ~ 100 kpc, we obtain a temperature of the outer, non- 
cooling region, T' to be 

T' ~ 1.3T, (21) 

for the typical values of Tc = 50 kpc, f3 = 0.6, r 25 oo = 600 kpc and rmax = 1 Mpc. If the 
temperature starts to decrease again at a larger radius, r ~ O.lSrigo as observed in nearby 
relaxed clusters (e.g. Vikhlinin et al. 2005; Markevitch et al. 1998), it would ahect the above 
equation by only < 10% due to a rapid decrease of the surface brightness at such large 
radii. Therefore T is likely to underestimate the ambient temperature by about 30% for 
Aooi ^ tage- Notice that this estimation is also consistent with the diherence of 30% in the 
average temperature between two subgroups shown in the previous subsection. Thus we 
consider that it is worth examining the luminosity - “ambient temperature” relation. 

By simply adopting the correction T' = 1.3T for 26 clusters with log(fcool/tage) < —0.5 
otherwise T' = T for 43 clusters, we plot the Ly^ — T' relation in Fig. 10c. As listed in Table 1, 
the difference in the normalization factors between r ^ 100 kpc mostly disappeared in this 
case. The slope of the power-law function is 2.80j]g;24 and thus marginally less steep compared 
with that of Lx—T. Fukazawa et al. (2004) recently derived scaling relations using the ASCA 
archival data on ~ 300 objects over a wide range of mass scales from elliptical galaxies to 
clusters. They found a smaller index of the Lx — T relation, a = 2.34 ± 0.29 at kT > 4 keV. 
Their analysis is different from the present one in that they utilized an integration area 
extending to 2-3 times rsoo for luminosity and the area excluding the central 1' — 3' region 
for temperature. Their temperature should therefore be closer to T' than T. In fact, the 
index of the Lx — T' relation mentioned above is slightly larger than but consistent within 
error bars with the result of Fukazawa et al. (2004). 

In Fig. lOd we show Likev(= Lx/T'"^'^^) versus (AcoiAage)- To quantify the dispersions 
and demonstrate their changes, we plot histograms of Tikev for the two Tc groups and those 
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with the above correction in Fig. 12. Since we are interested in the statistical signihcance of 
the change of the width, we performed the Gaussian htting with the method to include the 
Poisson error of each bin. The resultant Gaussian parameters (the mean, fi and the width, 
a) and the la errors are listed in Table 2 and also shown in each panel of Fig. 12. We thus 
found that by correcting the temperature with Eq. 21, the width signihcantly decreased from 
a/fi= (7.0T1.2) X10“^ to (3.3±0.9) x 10“^ (la error) for the smaller Tc group (see Fig. 12a, c) 
and the difference of the widths between two Vc groups disappeared: a/jj, = (3.3 ±0.9) x 10“^ 
and (3.5 ± 0.5) x 10“^ (Fig- 12c, d). Note that the la error of a//i is dominated by that of 
a and thus calculated by simply propagating the errors of a and fi. 

It is also remarkable that after the correction, there is no signihcant difference in the 
dispersions of the data points between fcooi ^ Gge and Tikev can be regarded as constant 
within the data scatter over two orders of magnitudes of fcooi (see also Table 2). Thus we 
consider that this also gives an evidence against the standard cooling-flow model because 
it predicts cooling that is increasingly enhanced with time. Our result indicates that the 
luminosity, that is the rate of thermal energy loss by the radiative cooling, does not depend 
on time very much, and some steady-state is attained for the gas of many small-core clusters. 
We will discuss about this in more detail in § 7. 

6. Lx - T/9 RELATION AND X-RAY MORPHOLOGY 
6.1. Lx — T/5 Relation 

In the case that the intracluster gas follows the isothermal /5-model, the virial temper¬ 
ature should be represented by Tvir = PT^ssX^il + x^)~^ ~ /5Lgas, where x = rvir/rc and 
typically a; 3> 1 (Ota 2001; Akahori & Masai 2005). We then replace T with /5T (here (3 is 
the slope parameter determined from the ROSAT HRI radial prohle htting) and derive the 
Lx — Tj3 correlation (Fig. 13a). 

As a result, the power-law index of the Lx — Tj3 relation is determined to be 2.67 ± 0.44 
, which is marginally smaller in comparison to that of the Lx — T. As for the data scatter 
around the mean relation, it becomes smaller for the small-core clusters: the width changes 
from a//i = (7.0±1.2) x 10“^ to (4.6±0.7) x 10“^ from Lx —T to Lx —T/5, while it becomes 
marginally larger for the large-core clusters: (4.2 ± 0.7) x 10“^ to (10.8 ± 7.8) x 10“^. The 
results are also summarized in Tables 1 and 2. 

To look into the cause, we investigated parameter correlations, (3 — T and f3 — Tc- We 

found that the (3 — T correlation is not signihcant (the correlation coefficient is 0.24) and 
-- 1 - 0.08 

at most /5 oc T -0.09. On the other hand, there is a trend to hud larger Vc for larger [3 
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noticeably at Vc > 100 kpc though such a trend is not clear for Vc < 100 kpc. In addition, 
several clusters show highly inhomogeneous surface brightness distributions and very large 
cores of > 400 kpc, which form a small peak at around ~ 500 kpc in the histogram 
(Fig. 5b). Considering that a large fraction of large-core clusters (24 out of 34) exhibit 
irregular X-ray surface brightness distribution, we may say that the above change in the 
dispersion likely reflects the X-ray morphology (see § 3.2 and Table 2 of Ota & Mitsuda 
(2004) for classihcation of the X-ray morphology). 

Particularly, all the hve clusters that are located far below the mean relation (i.e. ^30 
A1895, #50 AC118, #18 MS0451.5+0250, #46 1E0657-56, and #44 A483) have multiple 
peaks or a signature of cold front in the X-ray images, which can be attributed to cluster 
merging. In addition, they have very large core radii of > 400 kpc and large (3^1 
in comparison to the average value for the distant sample, (P) = 0.64 (Ota & Mitsuda 
2004). In Fig. 13b, they appear at the largest end of Tc- If we exclude the five clusters and 
estimate the width of the distribution of Tikev(= Lx/ikTpp-^'^), it results in a smaller value 
oi a/jj, = (6.4 ± 5.5) x 10“^ for 29 large core clusters though it is associated with the large 
statistical uncertainty. To show the relevance to the X-ray morphology in an alternative 
way, we divide the sample into the regular and irregular clusters and plot them by using 
different colors in Figs. 13c and d. We roughly see three phases: i) irregular clusters with 
very large core radius and smaller Likev at logrc > —0.4, ii) a coexistence of regular and 
irregular clusters with large core radius at —1 < logrc < —0.4, and iii) regular small-core 
clusters with higher Tikev at logr^ < —1. Note that for iii), there is only one exception, 
which is #17 A115, an irregular cluster lying at the smallest Tc- ^ Therefore we suggest that 
the above results show a trend that merging clusters particularly in the phase i) segregate 
from more relaxed clusters on the Lx — TP plane. 


6.2. Lx — T'P Relation 

In Fig. 14, we hnally show a plot of Lx — T'P relation for T' corrected in the same 
manner as in § 5.2 and Fig. 13. The results of the power-law fitting and the dispersions 


^A115 (z = 0.1971) exhibits two peaks in the surface brightness distribution and the temperature variation 
measured with ASCA implied that the cluster is indeed a merger system (Shibata et al. 1999). Gutierrez 
& Krawczynski (2005) reported based on the Chandra observation that the northern subcluster is highly 
nonuniform and there is a bright central core of about 10^^ in diameter, which might be explained by the 
ISM of the cD galaxy or a cooling core sloshing in the cluster. Our #model analysis with ROSAT yielded 
Tc = consistently with the Chandra result, and thus is very likely to have picked up the dense bright 

core of the northern subcluster. 
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around the mean relation are also listed in Tables 1 and 2. The power-law index is obtained 
to be 2.54lg;2g from the htting to 69 clusters. This is smaller in comparison with the Ly^ — T 
and other three relations derived above, yet signihcantly steeper than the self-similar model. 
In Fig. 14b, we hnd that the distribution of the data points is fairly constant against the 
cooling time: in particular for 43 clusters with fcooi < tage, the scatter is small: a/jj,= 
(2.7 ± 0.5) X 10“^ (Table 2). As suggested in the previous subsection, we see again the 
marginal trend of segregating the merging cluster from the regular clusters in Fig. 14a. 

In the next step, we further divide the sample into three subgroups according to the X- 
ray morphology: the regular single-/? cluster, the irregular single-/? cluster, and the regular 
double-/? cluster, and distinguish them by using different colors in Fig. 14c and d. From 
Fig. 14d, we may recognize a trend of the morphological change along the fcooi axis. Here 
it is possible to define three phases as follows: from larger to smaller tcooi, i) the irregular 
large-core clusters at log(tcooi/tage) ^ 0, where a large scatter of Likev is seen, ii) the regular 
large-core clusters at —0.3 < log(tcooi/fage) < 0, and hi) the regular small-core clusters at 
log(tcooi/tage) < “0.3. The regular and irregular clusters are no more mixed up in the 
.hikev — Aooi plot apart from a few exceptions (they are irregular clusters with small Vc, 
^17 A115, MS1910-I-6736, and ^11 MS0906.5-1-1110). It also should be noted that the 

boundary at logtcooi/fage = 0 in Fig. 14d, i.e. the boundary between i) and ii), approximately 
corresponds to the crossover point at logfcooi ~ login and logrc ~ —0.4 in Fig. 9a, and thus 
also to the boundary between i) and ii) in Fig. 13d. Unlike the Likev — ''"c plot (Fig. 13d) 
where the regular and irregular clusters are mixed in the intermediate range of Tc, the three 
groups are now clearly separated from each other in Fig. 14d. 

Furthermore, the double-/? clusters are located in iii) log(icooi/iage) ^ —0.3, the same 
region as the regular small-core clusters. As the inner-/outer-core dominant double-/? clusters 
are at a shorter/longer side of the range, they are suggested to be closely related to the small- 
/large-core regular clusters. We thus consider that this may be connected with the evolution 
of the X-ray surface brightness distribution after the onset of radiative cooling. We will 
discuss on this point in more detail in the next section. 


7. DISCUSSION 
7.1. Summary of the Results 

From the statistical analysis of 69 distant clusters in 0.1 < < 0.56 with the ROSAT and 

ASCA X-ray catalog (Ota & Mitsuda 2004), we have obtained the following observational 
results. 
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1. Regarding the cluster outer region, the density prohles are nearly universal and the 
scatter around the Lx{> 0.2r5oo) —T relation is surprisingly small. However, the slope 
of the relation is signihcantly larger than the self-similar expectation of 2. 

2. For the central region, the gas density exhibits a signihcant scatter and the self-similar 
condition is not satished particularly in the small core clusters. 

3. We studied the redshift evolution of the Lx — T relation. If restricted the luminosity 
integration area to the central 0 . 2 r 5 oo region, there might be a weak, negative evolution, 
however, the normalization factor of the relation does not signihcantly change with 
redshift and is consistent with the nearby sample within the data scatter. 

4. We investigated the parameter correlations focusing on the Lx — T relation and the 
connections to and tcooi, and suggested based on the X-ray fundamental plane anal¬ 
ysis that tcooi is likely to be a control parameter for the ICM structure of the cluster 
core region. 

5. For all small core clusters, tcooi < 4ge and the central electron density is systematically 
higher as Ueo oc which is suggestive of a central concentration of the gas as cooling 
progresses. On the other hand, we conhrmed that the temperature decrease is mild in 
the small core clusters; it is estimated as T = 1.3T' using the universal temperature 
prohle. Then we derived the luminosity-‘ambient temperature’ relation and showed 
that the dispersion of Li^ev signihcantly decreases and can be regarded as constant 
against tcooi within the data scatter. 

6. From the Lx — TjS and Lx — T'P relations, we showed that Likev is related to the X-ray 
morphology, and that there is a marginal trend for the merging clusters segregating 
from the regular clusters on their planes. We indicated the possible three phases of the 
gas property along the tcooi axis in Fig. 14d, according to the different level of cooling 
and ICM density structure. 

7.2. Impact of Radiative Cooling on the Lx — T Relation and Comparison to 

the Quasi-hydrostatic Model 

In the context of isobaric cooling how model, the local density of the gas increases 
against decreasing the temperature so that the thermal pressure P(r) at given r is kept 
constant within the cooling radius, i.e. P(r) = ne{r)kT{r) = const., and then Ue oc T~^ 
would follow. If this is the case, since Lx oc for the gas hotter than kT ~ 2 keV, the 

Lx — T relation would exhibit Lx oc For a typical cooling radius of rgooi ~ O.lrsoo 
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(~ 100 kpc), Lx(< Tcooi) contributes about 50% of the total luminosity. When we derived 
the Lx{< O.lrsoo) — T relation, however, we did not hud any tendency for such an inverse 
correlation (see also the Lx{< 0.2r5oo) — T relation in Fig. 4a). 

The observational results indicate that the gas in the core region of small core clusters 
is undergoing radiative cooling but the temperature decline towards the center seems mild. 
The gas may be considered close to a steady state. We will discuss such a possible state 
of the gas, quasi-hydrostatic cooling, proposed by Masai & Kitayama (2004). The model 
predicts a moderate and smooth gas inflow with hydrostatic balancing. In the context, unlike 
isobaric cooling flows that increase the local density so P{r) is kept against local cooling, 
quasi-hydrostatic cooling allows the gas to modify its prohle or core size so VP(r) matches 
with the force by gravitational potential. 

In § 5.2 we derived the Lx — T' relation by correcting the emission-weighted temperature 
for small-core clusters with short tcooi based on the idea of the universal temperature prohle, 
and showed that the bolometric luminosity normalized by the temperature dependence of 
Likev is constant within the data scatter over a wide range of tcooi/^age- This strongly 
suggests that the rate of thermal energy loss is kept nearly constant after the onset of cooling 
in the cluster core. 

According to the quasi-hydrostatic model, the temperature starts to decrease at the 
cooling radius and approaches a constant of ~ 1/3 the ambient temperature towards the 
center. The temperature prohle does not depend on its absolute value. This picture may 
explain the observations of nearby cooling-how clusters. The mass inhow rate from the 
outer region is controlled by cooling so as to maintain the quasi-hydrostatic balance and 
is expected to vary little through the how. The condition is well satished in the case that 
the temperature is kT > 2 keV where bremsstrahlung dominates cooling more than line 
emission. The present sample meets this condition since it has a temperature range of 
2 < kT [keV] < 12. Therefore our results can be consistently understood within a framework 
of the quasi-hydrostatic cooling model. 

At the same time, we have found that the gas density prohle for the outer part of the 
cluster does not vary very much from cluster to cluster, and is almost independent of the gas 
cooling at the center. This would also imply that the energy loss due to radiative cooling at 
the center is only a portion of the total energy of the thermal gas, and the outer region is 
likely to oher a great reservoir oi the heat, which is carried a fraction into the inner cooling 
region advectively by the gas inhow under quasi-hydrostatic balance. 

As described above, cooling with quasi-hydrostatic balancing modihes the core structure 
of the gas. This may give a clue to the origin of the small core of ~ 50 kpc or the inner core 
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component of the double-/? cluster. In addition the study would also have an impact on the 
present understanding of the origin of the dispersions around the Lx — T relation because 
it provided the observational evidence that they are directly linked to the thermal evolution 
of the ICM and showed how the relation would be modihed when the effect of the cooling is 
considered. 


7.3. Evolution of Density Structure of X-ray Clusters 

In this subsection we attempt to present an overall picture of the ICM evolution on the 
basis of the observational results. First we consider about the Tc distribution and the origin 
of the double-/? nature of the clusters from the phenomenological point of views. 

We have reported that the cluster core radius has a large cluster-to-cluster dispersion of 
more than an order of magnitude, even if restricted the sample to the regular clusters (Ota 
2001). This is partly due to the presence of two peaks separated by a factor of ~ 4. For 9 (7) 
out of the 35 regular clusters, the surface brightness is signihcantly (marginally) better htted 
with the double-/? model, and the average values of two core radii for the double-/? clusters 
coincide well with the two Tc peaks of the single-/? clusters. Then a question may arise: 
are there three subclasses of clusters, small-core single-/?, large core single-/?, and double-/? 
clusters, or are we looking at some different aspects of a single evolutionary process? From 
the imaging analysis, we noticed that the double-/? clusters are further subdivided into inner- 
core dominant and outer-core dominant ones; if we are to perform a single /? model htting 
for these clusters, we would pick up either of the core radii as the best-£t result. Thus the 
similarity in the distributions for the single-/? and double-/? clusters can be interpreted in 
the following way: the (regular) clusters have more or less double-/? nature in their structure, 
and when X-ray emission from one of the two cores dominates the whole cluster emission, 
we recognize it as the single-/? case and pick up either the inner or outer core as the core 
radius. This will account for an apparent large dispersion of core radius. Accordingly the 
single-/? clusters with a small/large core radius may respectively correspond to extreme cases 
of inner/outer core dominant double-/? clusters. 

Then why are there inner-core dominant and outer-core dominant clusters? Such double 
structures have often been attributed to the cooling-flow phenomenon at the cluster center. 
However, the present analysis showed that the standard cooling-flow model with dynamic 
inflow is not preferable to account for the results, in particular the T—Ueo and Likev—Aooi/^age 
relations. On the other hand, in § 6 we suggested the possible three phases of the ICM 
property on the Likev—Looi/^age plane in the light of the radiative cooling and the relevance to 
a variety of density structure. It is likely that the double-/? clusters are those in the transient 
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phase of their cooling cores, i.e., core size transition from outer-core dominated to inner-core 
dominated, mentioned in the quasi-hydrostatic cooling model. As the interpretation of the 
double core nature, we consider in the following way: when some time passed after the 
collapse and the cooling becomes important, the core is radiatively cooled and the central 
density gradually becomes higher to evolve from an outer-core-dominant cluster toward an 
inner-core-dominant cluster. 

The above discussion on the evolution of gas structure focused on the effect of radiative 
cooling and the possibility of quasi-hydrostatic state since the X-ray analysis showed that a 
majority of small core clusters are distributed in a narrow range of Tikev once the temperature 
bias is corrected (Fig. 12) and thus supports the existence of a typical value. Moreover, a 
long tail at the smaller Tikev is seen in the histogram for Tc < 100 kpc (Fig. 12c). This 
indicates existence of some peculiar clusters or extra physical mechanism that produces such 
a variation. Recently, O’hara et al. (2005) examined a relation between substructure and 
cool cores and the scatter about cluster scaling relations using 45 nearby clusters. They 
found that cool core clusters, which usually have less morphological substructure, exhibit 
higher intrinsic scatter about scaling relations even after correcting for cooling core effect, 
suggesting a possibility that a more global process is at work. Then the variation in the Ti^ev 
histogram for the small clusters may be related with the scatter found by O’hara et al. 
(2005). To further clarify both the structural evolution and the scatter of scaling relations, 
a systematic study with higher sensitivities will be important. 


7.4. Related Issues and Future Observations 

We will remark on the effect of the cold front detected in nearby clusters with the Chan¬ 
dra observations in regard to the present interpretation, and the assumption of the universal 
temperature prohle with consideration on the recent progress of Chandra and XMM-Newton 
observations. 

The cold front structure is characterized by a compact core with the density and tem¬ 
perature jumps in contrast to the surroundings, and interpreted as a remnant of subcluster 
merger (Markevitch et al. 2000). For higher redshifts, it is more difficult to directly identify 
such small structures due to the limitation of the sensitivity to measure the temperature 
maps, and some of them might possibly be fitted with the double-/? model. However, as 
the cold front is a signature of recent merger, it is expected to be accompanied with the 
irregular X-ray morphology. Thus for relatively bright clusters, a small core of the cold front 
in an irregular cluster may be distinguishable from that of a regular cluster by referring 
to the morphological information. For example, A115 and 1E0657-56, for which a possible 



24 


remnant of cooling core of a merging subcluster was pointed out (Gutierrez & Krawczynski 
2005; Markevitch et al. 2002), were classified into irregular clusters from our image analysis, 
though they have different Tc and are located in the different phases in Fig. 14d ( i) and iii), 
respectively). For more accurate identihcation of the cold front in distant clusters, it is nec¬ 
essary to measure the 2 dimensional density and temperature maps with higher sensitivity 
observations. 

We have utilized the universal temperature profile to estimate the effect of tempera¬ 
ture decrease and relate the emission-weighted temperature to the ambient temperature for 
clusters with short cooling time and apparent luminosity enhancement in Fig. 10b. Due 
to the limited spatial resolution of ASCA, the temperature prohles were not directly con¬ 
strained for the distant sample. However for some of them, the results have been already 
published with the Chandra and XMM-Newton data. As for the double-/? clusters or the 
small core cluster with very short cooling time of —1.3 < log (Gooi/Gge) ^ —0.7, A2204, 
PKS0745-191, A2390, A1835, the gradual temperature decline within the central region of 
roughly r < 100 — 200 kpc were reported (Sanders et al. 2005; Hicks et al. 2002; Allen et 
al. 2001b; Majerowicz et al. 2002). For the other double-/? cluster, A1689, with slightly 
longer logfcooi/Gge = —0.56, the nearly isothermal profile was found with Chandra (Xue & 
Wu 2002), and a complex temperature structure was further reported with XMM-Newton 
(Andersson & Madejski 2004). Therefore, from the overall standpoint, the assumption of the 
universal temperature prohle for logfcooi/Gge ^ —0.5 seems reasonable for many of the small 
core and double-/? clusters. To directly conhrm the temperature structure in relation to 
the gas cooling, we need to analyze the temperature prohle of individual clusters in various 
phases with using higher resolutions. 

In § 3.3, we derived the Tx(> 0.2r5oo)—T relation without correcting for the temperature 
prohle. Adopting the corrected temperature T' = 1.3T (Eq. 21), we obtain a similar slope, 
Lx(> 0.2r5oo) = 2.34+2-35 x for 69 clusters (y^/d.o.f. = 610/67) (Table 1) 

and a smaller standard deviation for Tikev, = (4-5 ± 0.5) x 10“^. Once the ambient 
temperature T' is directly measured for a large sample of clusters by future observations, we 
would suggest to use it instead of T in the analysis of cluster scaling relations. 
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A. A SIMPLE Lx-T-r-c MODEL 


In order to evaluate the Lx — T relation quantitatively, we assume the (3 prohle for the 
gas density; 

21 


eaSjO 


, r 

1 + - 


-3/3/2 


Then the total gas mass, Mgas, within r^ir is evaluated by 

Mgas = 47rr^pgas,o [ x'^{l + 


(Al) 


(A2) 


where x^ = r\i^/rc. For the observed range of a/v, i.e. 4 < a/v < 64 (see Fig. 1), the integral 
in Eq. A2 can be well approximated by a single power law function, thus is oc a;“ for a hxed 
value of (3 (Evrard & Henry 1991). Under the assumptions that the cluster mean density is 
Ac times the critical density of the universe, p = 3Mvir/47rr4j. = AcPcrit, and the gas-mass 
fraction, /gas = Mgas/Mvir is constant among clusters, Eq. A2 yields 

47r 

Afgas — AcPcrit/gas- (A3) 

Eqs. A2 and A3 relate the central gas density, Pgas.o, with Tc and as follows. 


Pgas,0 2 ^ 


3—a 


^cpcritf^ 


gas- 


The X-ray luminosity. Lx, 
Lx = 




is expressed as 

A(Tg„)4^r= ( 
A(Tg„)4irr/ ( 
1.435 X 10“^^ 

■/gas(AcPcrit) 


f x^{l-\-x^) ^^dx, 

m J Jo 

dgas,0 \ ^ 

rh J 

47r / 9(3k \ 

9m^ \47rmG' / 

■l-\-a—b/2rTi7/2—a+b/2 —3+2a—b 
^ gas ' c ’ 


(A4) 

(A5) 

(A6) 

(A7) 

(AS) 


where A is the emissivity of the free-free emission and rh = prrip {p : the mean plasma 
weight, 0.63). In the case of P = 0.7, a = 0.96, and b = 0.025 for x^ = 32 {x^ dependence is 
small), which leads 

Lx oc /gis(A.Pcnt)-°-°^Tg2a^e-^A (A9) 

We dehne the virial temperature as 


kT ■ — 

vir — 


mGMvir ATTUlGAcPcvitr'^ir 


3r, 
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(AlO) 
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and assume that it is related to the gas temperature through kT^i-c = P'kT^as- 

In Fig. 7, we show the Lx — T relation numerically calculated according to Eq. A8 for 
r'c = 50,100,200 /ifg^kpc, f3 = 0.7, and = 500. The gas-mass fraction within rsoo) 
is hxed at the mean value of /gas = 0.27 for the distant sample (Ota & Mitsuda 2004). We 
should notice that dependence of Lx — T relation on pcntAc is very small; the power-law 
index is —0.05, and is almost independent of the cluster-collapse epoch. Then the collapse 
redshift is assumed to be Zcoi = 0.3, which is an average observed redshfit of the sample. We 
adjusted the value of (3' to jS' = 0.4/3 in the calculation so as to make the absolute values of 
Lx consistent with the observations. 
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Fig. 1.— (a) Comparison of two luminosity measurements. The bolometric luminosities 
within rsoo estimated from the iJOFTT HRI are plotted against those from the ASCA CIS. 
The clusters are divided into three subgroups according to the best-£t temperature from the 
ASCA spectroscopy, 18 clusters with 2 < kT [keV] < 5, 30 clusters with 5 < kT [keV] < 7 
and 21 clusters with 7 < kT [keV] < 12 and denoted with the crosses, the asterisks, and the 
open circles, respectively. The solid line corresponds to Tx,hri = -f^x.cis- (b) Comparison of 
the cooling time, Rooi estimated with the ROSAT and ASCA analysis (Ota & Mitsuda 2004) 
and Chandra (Bauer et ah 2005) regarding 21 clusters with 0.15 < 2 ; < 0.35. 



0.01 0.1 1 


Fig. 2.— rsoo — Tc relation (a) and histogram of (b) for 69 clusters. In the panel (a), 35 
clusters with small core of Tc < 100 kpc and 34 clusters with larger core of Tc > 100 kpc 
are shown with the green asterisks and the magenta circles, respectively. The dotted lines 
indicate the self-similar condition corresponding to four different constant values of tsoo/tc- 






30 



radius [r^og] 


Fig. 3.— Electron density profiles for 69 clusters. The best-fit density profiles derived with 
the single /d-model are plotted, where the radius is normalized with O.lrsoo is indicated 
with the vertical dotted line, inside which the scatter is the most prominent. 



log(kT [keV]) log(kT [keV]) 


Fig. 4.— Tx(< 0.2r5oo) — T relation (a) and Lx(> 0.2r5oo) — T relation (b). The sample is 
divided into two Tc groups and thus the meaning of the symbols are the same as Fig. 2. The 
solid line represent the best-£t power-law for 69 clusters in each panel. The best-£t models 
obtained for Tc < 100 kpc and > 100 kpc are also shown with the dashed and dot-dash lines, 
respectively. 
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Fig. 5.— Lx — T relations for three subsets of the data: 18 clusters with 0.1 < z < 0.2 (blue 
crosses), 27 clusters with 0.2 < 2 : < 0.3 (green asterisks), and 24 clusters with 0.3 < 2 < 0.56 (red 
circles). In the panels (a) and (b), the luminosity was calculated within the overdensity radius rsoo- 
The solid line shows the best-fit power-law model for 69 distant clusters. In the panels (c) and (d), 
the Lx(< 0.2r5oo) — T and Lx(> 0.2r5oo) — T relations are shown, respectively. In the panels (e) 
and (f), the luminosity was calculated for a fixed radius of 2 Mpc. For 26 clusters with a short 
cooling time, we corrected the effect of possible cool emission regarding both the luminosity and 
temperature (see text for details). The solid line indicates the local relation obtained by Markevitch 
(1998). In the panel (f), the luminosity was further scaled with redshift as Lx(l -|- In each 

panel, the blue, green, and red dashed lines show the best-fit power-laws for the above three redshift 
ranges, respectively. In the panel (a), the slope of the relation was included as a free parameter 
in the fits, while in the panels (b)-(f), the slope was fixed at a = 3.01, 4.37, 2.85, 2.64, and 2.64, 
respectively. 
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Fig. 6. — Lx — T relation (a) and Likev — ''"c relation (b). In the panel (a), we find the 
smaller core clusters are brighter on average compared to the large core ones of the same 
temperature. In the panel (b), Likev(= Lx /is the normalization factor of the Lx—T 
relation derived for each cluster. The solid line indicates the best-fit normalization factor of 
the Lx — T relation for 69 clusters in the panel (b). The negative correlation between Tikev 
and Tc corresponds to the luminosity difference between the two subgroups, which is roughly 
approximated as Ti^ev oc r~^ (the dotted line). 


ITD 



Fig. 7.— Lx — T relations expected for clusters with a constant core radius, Vc = 50 or 100 
or 200 kpc (dotted lines). /? = 0.7 is assumed in the model calculation. See Appendix A for 
details of the model. 
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Fig. 8.— Results of X-ray fundamental plane analysis. In panel (a), the best-fit log X —log Z 
plane obtained for the distant clusters is shown. The meaning of the symbols are the same 
as Fig. 2. In the panel (b), 45 nearby clusters taken from Mohr, Mathiesen & Evrard (1999) 
were projected onto the same fundamental plane as panel (a). According to Table 2 of Mohr, 
Mathiesen & Evrard (1999) non cooling-flow clusters are shown with the red triangles, and 
inner-core and outer-core components of cooling-flow clusters are separately shown with the 
solid blue boxes and open black boxes. 
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Fig. 9.— fcooi — 'f'c (a) and T — Ueo (b) relations. For panel (a), tcooi is the radiative cooling 
timescale (see text). There is a strong correlation with Tc which follows fcooi oc The 
Hubble time, tu = 13.4 Gyr (Spergel et al. 2003) is indicated with the horizontal dotted line. 
Then tcooi < tn for all clusters belonging to the smaller core group. For panel (b), Ueo is the 
central electron density obtained from the (3 model analysis. There is no clear difference in 
the temperature range between the two subgroups though Ueo is scattered nearly over two 
orders of magnitudes. Their average temperatures of 5.4 keV and 7.0 keV are shown with 
the dashed line and dot-dash line, respectively. 
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log(kT’ [keV]) ^°S(*-cool/^age) 

Fig. 10.— Lx — T relation and tcooi- Upper panels, Lx — T relation of clusters (a) (the same 
as Fig. 6a), and Fikev as a function of fcooi/tage (b). Bottom panels, the luminosity to the 
‘ambient temperature’ relation, Lx—T' (c) and Likev = Lx/{kT'Y'^^ as a function of tcooi/tage 
(d). In the panels (c) and (d), for clusters with short cooling timescale of logfcooi/^ge < —0.5, 
the temperature decrease is corrected with T' = 1.3T, otherwise T' = T (see § 5.2). We hud 
a decrease of the dispersion around the mean Lx — T' relation in comparison to the Lx — T, 
which we show in a more quantitative manner in Fig. 12. 



l°g(*-ccol/tage) l°g(hcal/hge) 

Fig. 11.— Fikev as a function of tcooi Aage in three redshift ranges. In the panels (a) and (b), 
Uikev = Lx/ikTY''^^ (§ 5.2) and Lx/{kT'(3)“^-^^ (§ 6.2) are plotted, respectively. 
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Fig. 12.— Histograms of Fikev- The distributions of Likev(= L^l for 35 small core 
clusters and 34 large core clusters are shown with the green and magenta stepped lines in 
the panels (a) and (b), respectively. The results of Likev(= Ly^/{kT'Y'^^) are shown in the 
panels (c) and (d). The best-£t Gaussian model is shown with the dotted cnrve in each 
panel. The Poisson error of each bin was considered in the fitting thongh it is not displayed 
in the panels. The fitting results are summarized in Table 2. 
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Fig. 13.— Lx — T(3 relation and X-ray morphology. Upper panels, Lx — T(3 relation of 
clusters (a) and Likev = Lx/{kT as a function of fcooi/tage (b). Bottom panels, the 
same as panel (a) (c), and Likev as a function of Tc (d). The meaning of the symbols for two 
Tc classes are the same as (a) and (b), but the regular and irregular clusters are shown in 
blue and red respectively. In the panel (d), we indicated the three phases corresponding to 
i) the irregular clusters with a very large and a signature of merger, ii) the coexistence of 
the regular and irregular clusters with a large Tc, and hi) the regular clusters with a small 
Tc (see § 6.1). 
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Fig. 14.— Lx — T'j3 relation and X-ray morphology. The double-/? clusters are denoted 
with the light-blue stars. The hlled (open) stars correspond to the inner-core (outer-core) 
dominant double-/? clusters. In the panel (d), we indicated the three phases according to the 
different level of cooling and density structure: i) the large core, irregular clusters with long 
tcooi relative to tage, h) the large core, regular clusters with moderate tcooi, and iii) the small 
core, regular clusters with short tcooi- In the phase iii), 9 of the sample shows the significant 
double-/? structure, and the inner-core/outer-core dominant double-/? clusters are located at 
the shorter/longer cooling time (see § 6.2). 
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Table 1: X-ray luminosity-temperature relations for the distant sample. 




Sample 

Relation“ 

x2/d.0.f.’' 

Lx-T 


all (69) 

7-4ll4;40 X 

1603/67 



rc < 0.1 Mpc (35) 

1.05l(J;ii X 10‘‘3(^y)3.oi 

722/34 



Tc > 0.1 Mpc (34) 

4.52;°-37 X 

172/33 



rc < 0.1 Mpc (35) 

1.45;j;g X 1043(fcT)2-827S:f4 

716/33 



rc > 0.1 Mpc (34) 

4-37;2;3o X 1042(fcT)^ °^-o-3t 

172132 

Lx-r 


all (69) 

7Altai X 10-^2 

669/67 



rc < 0.1 Mpc (35) 

8.041°;^^ X l0'^\kT'^ 

461/34 



Tc > 0.1 Mpc (34) 

6-76;(Jj^ X 1042(^T')^ ®° 

177/33 



Tc < 0.1 Mpc (35) 

1.04t°;3| X 1043(^y/)2.67t°;|^ 

457/33 



rc > 0.1 Mpc (34) 

4-37t2;30 X 1042(^y/)3.03t°;|| 

171/32 

Lx-Tfi 


all (69) 

8.32+3;32 ^ io43(/jy^)2.67l“::: 

1310/67 



rc < 0.1 Mpc (35) 

1.08l°;?® X 104‘‘(A:Td)2-67 

420/34 



rc > 0.1 Mpc (34) 

3.35l(J;4^ X 10‘^3(^y^)2.67 

169/33 



rc < 0.1 Mpc (35) 

1.10l)J;4^ X 

420/33 



rc > 0.1 Mpc (34) 

2.85l^;^8 X 1043(^y^)2.77lgJ^ 

169/32 

Lx - ri3 


all (69) 

e.ioli;^^ X 10‘‘3(^y/^)2.54l“;- 

628/67 



rc < 0.1 Mpc (35) 

Qmtia X 10‘^3(^y/^)2.54 

294/34 



rc > 0.1 Mpc (34) 

4.03l[!i[; X 10‘‘3(^y/^)2.54 

172/33 



rc < 0.1 Mpc (35) 

7.161^29 X 1043(/cT'/3)2-3i7S:28 

293/33 



rc > 0.1 Mpc (34) 

2.85t^;?| X 10‘^3('^y/^)2.77t);:®^ 

169/32 

Lx{< 0.2r5oo) 

-T 

all (69) 

3.021^^8^ X 10-^1 (fcr)^-3^-i'-5o 

2995/67 



rc < 0.1 Mpc (35) 

5.82l(;?^ X 1041(A:T)4-37 

1289/34 



rc > 0.1 Mpc (34) 

8.32+X38 X 10‘^°(A:T)4-37 

301/33 



rc < 0.1 Mpc (35) 

2.88+2^3™ X 1042(fcr)3-43t;j:^| 

1231/33 



rc > 0.1 Mpc (34) 

9.23;6°7f X 104i(fcr)3-^^-o'"9 

252/32 

Lx{> 0.2r5oo) 

-T 

all (69) 

4.95^1^3 X 1042(A:T)2-33lo:27 

551/67 



rc < 0.1 Mpc (35) 

S-Selase X lO^^^^T)^-®^ 

321/34 



rc > 0.1 Mpc (34) 

4.27l°j® X 10‘^2(^y)2.85 

174/33 



rc < 0.1 Mpc (35) 

7.4i;^-^® X 1042('^y)2.68;«;« 

317/33 



rc > 0.1 Mpc (34) 

1.45;J;^5 ^ io42(fcy)3.4o;«:- 

148/32 

Lx{> 0.2r5oo) 

-T' 

all (69) 

X 10-^2(^y/)3.057“;99 

610/67 



rc < 0.1 Mpc (35) 

2.04+°;23 X 10-^2(^y/)3.05 

362/34 



rc > 0.1 Mpc (34) 

2.88l(J;26 X 10-^2('^y/)3.05 

158/33 



rc < 0.1 Mpc (35) 

5.07l^j? X 10-^2(^y/)2.59+°;« 

336/33 



rc > 0.1 Mpc (34) 

7A5tlil X 1042(fcT')^ '‘°-o'« 

148/33 


“The power-law relations derived for the distant sample. The results for two Vc subgroups with/without fixing 
the slope at the best-fit value obtained for 69 clusters are also shown. The measurement uncertainties of 
both the X and y axes are took into account in the fit. 
and the degree of freedom of the power-law fitting. 
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Table 2: Means and standard deviations of logLikev- 


log T IkeV 

Sample 



xVd.o.f.^ 

a//i [10-3]^ 

logLx/(A;T)3-oi 

all (69) 

42.75 ± 0.04 

0.31 ±0.04 

4.73/5 

7.3 ±0.9 


re < 0.1 Mpc (35) 

42.92 ±0.05 

0.30 ±0.05 

1.95/5 

7.0 ± 1.2 


Tc > 0.1 Mpc (34) 

42.60 ±0.04 

0.18 ±0.03 

0.82/2 

4.2 ±0.7 


^cool ^ ^age (43) 

42.88 ±0.05 

0.30 ±0.04 

2.07/5 

7.0 ±0.9 


^cool ^ ^age (26) 

42.57 ±0.04 

0.18 ±0.04 

0.56/2 

4.2 ±0.9 

logLx/(A;T')2-80 

all (69) 

42.82 ±0.03 

0.19 ±0.02 

4.53/5 

4.4 ±0.5 


re < 0.1 Mpc (35) 

42.92 ±0.03 

0.14 ±0.04 

7.35/5 

3.3 ±0.9 


re > 0.1 Mpc (34) 

42.78 ±0.03 

0.15 ±0.02 

0.78/2 

3.5 ±0.5 


^cool ^ ^age (43) 

42.87 ±0.04 

0.19 ±0.04 

5.23/5 

4.4 ±0.9 


^cool ^ ^age (26) 

42.77 ±0.04 

0.16 ±0.03 

0.99/2 

3.7 ±0.7 

logLx/(A;T/3)2-67 

all (69) 

43.63 ±0.05 

0.35 ±0.05 

11.53/6 

8.0 ± 1.1 


re < 0.1 Mpc (35) 

43.92 ±0.04 

0.20 ±0.02 

7.67/3 

4.6 ±0.7 


re > 0.1 Mpc (34) 

43.69 ±0.76 

0.47 ±0.34 

0.75/3 

10.8 ±7.8 


^cool ^ ^age (43) 

43.86 ±0.04 

0.22 ±0.03 

11.22/4 

5.0 ±0.7 


fcool ^ ^age (26) 

43.44 ±0.27 

0.41 ±0.20 

2.62/3 

9.4 ±4.6 

logLx/(A:T'/3)2-54 

all (69) 

43.60 ±0.27 

0.26 ±0.20 

18.98/4 

6.0 ±4.6 


re < 0.1 Mpc (35) 

43.83 ±0.02 

0.11 ±0.02 

3.58/3 

2.5 ±0.5 


re > 0.1 Mpc (34) 

43.81 ±0.02 

0.29 ±0.06 

2.66/3 

6.7 ± 1.4 


^cool ^ ^age (43) 

43.81 ±0.02 

0.12 ±0.02 

5.92/3 

2.7 ±0.5 


^cool ^ ^age (26) 

43.41 ±0.07 

0.29 ±0.07 

2.13/3 

6.7 ± 1.6 


“The Gaussian mean, /i and width, a obtained from the fitting to the TikeV distribution. The quoted 
errors are the Ict. 

^X^ and d.o.f. of the Gaussian fitting. 

“a/zi and the Ict error in 10“^. 



